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STARTUP ANALYSIS FOR A HIGH -TEMPERATURE GAS-LOADED HEAT PIPE 

by Peter M. Sockol 
Lewis Research Center 

SUMMARY 

A model for the rapid startup of a high -temperature gas-loaded heat pipe is pre- 
sented. A two-dimensional diffusion analysis is used to determine the rate of energy 
transport by the vapor between the hot and cold zones of the pipe. The vapor transport 
rate is then incorporated in a simple thermal model of the startup of a radiation -cooled 
heat pipe. Numerical results for an argon-lithium system show that radial diffusion to 
the cold wall can produce large vapor flow rates during a rapid startup. The results also 
show that startup is not initiated imtil the vapor pressure p^ in the hot zone reaches a 
precise value proportional to the initial gas pressure p^. Through proper choice of p., 
startup can be delayed until p.^ is large enough to support a heat -transfer rate sufficient 
to overcome a thermal load on the heat pipe. 


INTRODUCTION 

A gas -loaded heat pipe contains a fixed quantity of a noncondensible gas in addition 
to its working fluid. The introduction of a gas into an ordinary heat pipe has been sug- 
gested as a means of aiding its startup. When a high -temperature heat pipe is closely 
coupled to a thermal sink, the very low initial vapor pressure may make it incapable of' 
delivering the energy demanded during startup (ref. 1). Some technique must be found to 
decrease the coupling between the pipe and the sink if startup is to be successful. 

During operation of a gas -loaded heat pipe, vapor flows from the evaporator to the 
condenser and the noncondensible gas is carried with it. The gas collects in the end of 
the condenser and decreases the effective heat -rejection area. As the temperature, and 
hence the v£q)or pressure, in the evaporator is increased, the gas is compressed to a 
smaller volume. If a siifficient quantity of gas is introduced, the heat pipe can be com- 
pletely decoupled from the thermal sink until the vapor pressure is high enough to sup- 
port large energy transport rates. 



Before a gas -loaded heat pipe is incorporated in a larger system it is necessary to 
be able to predict its performance during startup. To this end a two-dimensional anal- 
ysis of the diffusion of vapor through the gas in the condenser is performed in the present 
work. The energy transport rate obtained in this study is combined with a simple ther- 
mal model of the startup of a radiation -cooled heat pipe. Numerical results are pre- 
sented for a long lithium heat pipe with argon as the noncondensible gas. They shovild be 
of value in the design of experiments with high temperature gas-loaded heat pipes. 
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SYMBOLS 

expansion coefficient, eq. (18) 
parameter, eq. (31) 
parameter, eq. (32) 
expansion coefficient, eq. (23) 
expansion coefficient, eq. (19) 
expansion coefficient, eq. (24) 
heat capacity per unit length 
total molar density 

initial molar concentration of noncondensible gas 

outside diameter of heat pipe 

parameter, eq. (36) 

parameter, eq. (37) 

binary diffusion coefficient 

exponential integral 

fimction, eq. (16) 

function, eq. (17) 

dimensionless vapor flux 

heat of vaporization per unit mass 

Bessel fimction 

length of hot zone 

length of evaporator 
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Pi 

Py(T) 
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Qv 
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Tc 
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length of heat pipe 
dimensionless L 
dimensionless 
dimensionless Lp 
molecular weight of vapor 
vapor flux, eq. (3) 
molar vapor flux 

initial pressure of noncondensible gas 

vapor pressure 
heat input 

heat lost by radiation 

energy transported by vapor 

dimensionless Q 

dimensionless 

dimensionless 

radius 

inside radius of wick 
dimensionless R 
temperature of hot zone 
temperature of cold zone 
time 

time at which Q is switched on 

mole fraction of noncondensible gas 

mole fraction of vapor 

concentration variable , eq. (9) 

distance along pipe 

dimensionless Z 

zeros of J ' 
o 

zeros of 
Euler constant 



6y Kronecker delta, eq. (20) 

e emissivity of pipe wall 

0 dimensionless T 

X(e,l) function, eq. (8) 

X parameter, eq. (33) 

U 

cr Stefan -Boltzmann constant 

Subscripts: 

c cold zone 

h hot zone 

Superscript: 

( ) ' differentiation with respect to argument 


DESCRIPTION OF THE MODEL AND ASSUMPTIONS 

! I 

Thermal Model ' 

The thermal model has already been applied to the startup of an ordinary radiation- 
cooled heat pipe (refs. 2 and 3). When the power is txirned on, a hot zone of uniform wall 
temperature T forms in the evaporator. Once the vapor pressime p^(T) is high enough, 
the front edge of this zone begins to move into the condenser. Despite the steep wall- 
temperature gradient at the front, it is assumed that axial heat conduction in the wall is 
negligible compared to energy transport by vapor diffusion. The remaining cold zone re- 
tains its initial temperature T . 

, This model is in definite opposition to the results of the steady -state analysis of ref- 
erence 4. There it is shown that at the front axial heat conduction in the wall is much 
greater than energy transport by the diffusing vapor. The present calculations confirm 
that the vapor flow rate past the front becomes quite small as steady state is approached. 
They also show, however, that very large diffusion rates accompany the motion of the 
temperature front during startup. Such rates are incompatible with a stationary front. 


Diffusion Model 

The goal of the diffusion analysis is to determine the vapor flow rate past the tern - 
perature front as a function of the length and temperature of the hot zone. To simplify 


4 



the calculations the local condensation rate upstream of the front is assumed to be neg- 
ligible, whereas the cold wall is taken to act as a perfect sink for the vapor. The transi- 
tion region between these two boxmdary conditions is considered to be vanishingly thin. 
The nonuniform vapor velocity in the hot zone must produce a secondary flow of noncon- 
densible gas. This effect is neglected, however; the upstream velocity is assumed uni- 
form over the cross section, and the gas is taken to be stagnant throughout the pipe. The 
axial pressure drop is neglected, and the total pressure, gas plus vapor, is set equal to 
the vapor pressure p^(T) in the evaporator. This restricts the analysis to cases where 
the vapor velocity is much less than the speed of sound during startup. In a long cylin- 
drical heat pipe the vapor concentration should relax to a quasi-steady profile in a time 
of order (R^/2. 4) /9, where is the inside radius of the wick and 9 is the binary 
diffusion coefficient for the vapor -gas mixture. If the distance moved by the temperature 
front during this time is small, then steady -flow equations can be used in the diffusion 
analysis. In an argon -lithium heat pipe at 1300 K, with = 1 centimeter, the relaxa- 
tion time is about 10 milliseconds. For a typical case (fig. 3) the maximum frontal mo- 
tion during this time is 0. 02 centimeter which is small. Finally, the product cQ, where 
c is the total molar density, is assumed constant and evaluated at the average of the 
temperatiu-es T and T^. 

DIFFUSION ANALYSIS 

t t 

The molar flux of vapor through a stagnant gas is given by (ref. 5) 

- XyN^ - V x^ = c^ V In Xg (1) 

where x^ and Xg are the mole fractions of vapor and gas, respectively. Conservation 
of mass requires 



In Xg = 0 (2) 

A cylindrical coordinate system is introduced with the main flow in the Z direction and 
the origin at the interface between the hot and cold zones (fig. 1). It is assumed that the 
axial distance over which Xg departs significantly from 0 or 1 is small compared to L 
or Lp - L where L is the length of the hot zone and Lp is the length of the heat pipe. 
Hence, the boundaries on equation (2) are taken at Z = ±oo. The boimdary conditions for 
equation (2) are as follows: 
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R = 
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1 

m ' ' 

^ X\\ 

Hot 

Cold 

Figure L - Coordinate system for diffusion analysis with main flow 


in Z-direction and origin at interface between hot and cold zones. 


For Z < 0: 


c ^ — In as Z — -<» 

az s o 


(3) 


-i In x„ = 0 at R = (4) 

3R S ^ 


For Z > 0: 


Xg — 1 as Z ~ oo 


(5) 


Xg = 1 at R = R^ 


( 6 ) 


where is an unknown constant flux. Conditions (4) and (6) represent no condensation 
and the perfect sink, respectively. The unknown is determined by specifying the ini- 
tial molar concentration Cj^ of the noncondensible gas at temperature T^. A mass bal- 
ance on the gas gives 


c.L 


P 



(7) 


where the bar indicates an average over the cross section, and Cj^ and c^ are constants 
obtained by evaluating c at the temperatures T and T^, respectively. Note that the 
integrands in equation (7), as finally written, go rapidly to zero away from Z = 0. 
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Dimensionless variables are introduced; 



( 11 ) 

( 12 ) 

( 13 ) 

( 14 ) 
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Equation (7) becomes 


F^(g) Fj^(g) =X(0,Z) 


with 


F^(g) = 2 



(1 - e“^)dz 


Fn(g) 


r r 

-21 r dr / 

Jn J-oo 


-gy 


dz 


(15) 


(16) 


(17) 


Once y(r,z) is determined, equation (15) can be solved numerically to obtain g, and 
hence the mass flux as a function of T and I . 

The solution of equations (10) to (14) is given by the following: 

z < 0: 

00 

y = Jo'(«i) = o ( 18 ) 

i=l 

z > 0: 

<XJ ^ 

y =22®i® Jo(%) = 0 (19) 

i=l 


Note that 0 !j = 0 and J^(0) = 1. The Bessel functions J^(o'.r) and J^(;3^r) satisfy the 
following relations (ref. 6, ch. 11): 
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where 6- = 1 and 6^ = 0 for i j. Requiring y and dy/dz to be continuous at z = 0 
and applying relations (18) to (20) to the results give 

00 

a. = -2^b.i3.(cy2 - j = 1 , . . . , oo (23) 

i=l 

OO 

b. = 2^a. ^ ^^-2^ j = 1, . . . , (24) 

i=l 

where a^ = AjJ^(or.) and b^ = -BjJj(/3|). A finite subset of equations (23) and (24) with 
j = 1, . . . , n can be solved by iteration. Final values are obtained by increasing n un- 
til the low order coefficients no longer change. For n = 30 the first few coefficients are 


aj = 0.2537 

bj = 0.2578 

a2 = -0.0761 

bg = 0.0638 

ag = -0.0336 

bg = 0.0308 

a^ = -0.0201 

b^ = 0.0188 

ag = -0.0138 

bg =0.0129 


The remaining coefficients decrease quite slowly, but the exponential factors in equa- 
tions (16) and (17) ensure adequate convergence away from z = 0. 

Once g has been found by solving equation (15), the energy transport by diffusion 
from the hot to the cold zone is given by 

Qv = V 

where is the molecular weight of the vapor and hj^ is its heat of vaporization per 
unit mass. For g of order one this corresponds to a very small energy transport. As 
an example, in an argon-lithium heat pipe at 1300 K with R = 1 centimeter and g = 1, 
equation (25) gives a; 7.4 watts. For large transport rates to occur it is therefore 
necessary to have g » 1. From equations (9) and (19) this corresponds to the bulk of 
the gas being pushed far enough into the cold zone to enable the vapor to reach the cold 
surface by diffusing radially through a very thin layer of gas. 

The functions Fg(g) and Fj^(g) must be evaluated numerically. The results have 
been approximated by analytical expressions for subsequent use. For g » 1, however, 
F^ is dominant and the following approximation is of interest. If only the first term in 
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the series of equation (19) is retained, equation (16) gives 


Fc(g) ^ ^ Iq -SBj^J^(/3^r)e |dz 


gBjJ^(^ir) 


(1-e-t)^ 


/'rdrsjl / 
■^0 ^ ^0 


2 y* r dr i3-l[y + In [gBjjQ(0jr)]+Ej[gBjJ^(^ir)]J 


,-l 


(In g - 1.339) 


(26) 


where y = 0. 5772. The radial integral was performed numerically. Note that the expo- 
nential integral (ref. 6, ch. 3) Ej(x) approaches zero for large x. Hence, F^(g) has only 
a weak dependence on g for g » 1. 


TIIVIE-DEPENDENT THERMAL ANALYSIS 

The time -dependent equations are written for a radiation-cooled heat pipe with a 
shielded evaporator. The state of the heat pipe at time t is shown in figure 2, where Q 
is the heat input to the evaporator, the heat lost by radiation, and the energy 
transported by the vapor to the cold zone . 

When the hot zone is expanding into the condenser, energy balances on control vol- 


4 . 
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umes enclosing the hot and cold zones, respectively, at time t give 


CLlii = Q-Qj. -Q^ (27) 

dt 

C(T-TJ^ = 0 (28) 

dt ^ 

where C is the heat capacity per unit length of the wall, wick, and liquid. In terms of 
0 and I these equations become 


Z ^ = q - q - q (29) 

dt ^ ^ 

(0 - 1) ^ = q (30) 

dt 

with q = Q/CR T . The quantities q_ and q are given by 

W L/ i V 

Qj. = Aj.0^ {I - Zg) (31) 

\ = V (32) 

where = irDeaT^/C , Zg = L^/R^, and = ;rM^(c^)hjg/CTg. Also, D is the outside 
diameter of the pipe , e the emissivity, and a the Stefan -Boltzmann constant. The for- 
mulation is completed by equation (15) with X(0,Z) given by equation (8). 

During the initial startup period, q^ remains negligibly small until Py(T) > p^^ and 
X > 0. For X > 0, there is an extremely rapid transient in which q^ approaches q 
and d0/dt becomes small. Thereafter, both Z and 0 increase at moderate rates, and 
q^ decreases steadily as q^. increases. This continues until steady state is approached. 
The behavior is illustrated in figure 3 for an argon-lithium heat pipe. 

In this latter period, equation (15) causes a severe numerical difficulty. The fimc- 
tion X(0,Z) has a very strong dependence on 0 and Z, while F^(g) is a weak function of 
g for g » 1. Thus, q^ is extremely sensitive to small changes in 0 and Z, which 
makes equations (29) and (30) numerically unstable for any reasonable size integration 
step (>0. 1 sec). 

Since Z(d0/dt) is small compared to (q - q^.) in this period, equation (29) can be used 
to obtain an approximate value of q^ and hence g. Equations (30) and (15) can then be 
solved simultaneously for Z and 0 . It is necessary, however, to first obtain a good ap- 
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Temperature. 









Time. t. sec 

(c) Initial gas pressure, pj, 200 torr; time at which heat input is turned on, t|, -13. 9 seconds. 

Figure 3. - Concluded. 

proximation to dfl/dt. Differentiating equations (15) and (8) gives 

K- F. (g)l = rF 'Cg) - 0-^ F'(g)l ^ (33) 

L ® " J dt dt L ^ J dt 

9 

with X. = T /p . Equation (29), with Z(de/dt) assumed small and q taken as 

C7 C IV V 

constant, gives 

- ^ {I -i)^ + ^1 (34) 

dt L ® dt dtJ 

Equations (33) and (34) give 

D- - D, ^ « 0 (35) 

® dt ^ dt 

with 

\ " (f'c - J'h)^ • - 'e' »6) 
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(37) 



Equations (29), (30), and (35) give 


Z + (0 - 1) 


D, 


— q - Oj 

dt ^ 


Finally, equations (29) and (38) give 


(38) 


(q - qr.) 


(0 - 1)D 


(0 


0 

1)D0 ^ 


(39) 


Equations (30), (15), and (39) form a system of one differential and two algebraic 
equations whose solution approximates that of equations (29), (30), and (15) when d 0 /dt 
is small. The stability of this approximate set of equations is much improved over that 
of the original set. 

Both sets of equations have been incorporated in a numerical integration scheme. 
The initial conditions at t = 0 are taken to be = Iq and X(0 q, Iq) =0. From equa- 
tion ( 8 ) this gives 


At earlier times q^. =0, q^ 0, and 0 = 1 + q(t - t^)/Zg where t^ is the time when Q 
is switched on. The integration is started using the original set of equations and a very 
small integration step. When q^ is close to its maximum value, the integration is con- 
tinued using the approximate equations. The numerical routine automatically varies the 
size of the integration step to achieve the maximum value consistent with accuracy and 

_4 

stability. In a typical case the stepsize varied from 10 second at the start to 10 sec- 
onds at the end of the integration. In a few cases the approximation was checked by using 
the original set of equations throughout the calculation. The agreement between the two 
methods was better than three significant figures in I and 0 . Solutions have also been 
obtained for a startup in which Q is increased in steps until the final steady state is 
reached. 
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NUMERICAL RESULTS 


The time -dependent theory of the preceding section has been applied to the startup of 
an argon -lithium heat pipe. The specifications of the pipe are as follows: Lp = 300 cen- 
timeters, L = 30 centimeters, D = 1.9 centimeters, R = 0.78 centimeter, C = 1. 99 
c w 

joules per centimeter per degrees Kelvin, and €=0.3. The values for R^, C, and e 
were obtained by selecting a 1. 0 millimeter thick tantalum alloy wall, an 0. 5 -millimeter - 
thick annular liquid return, and an 0. 2 -millimeter -thick tungsten screen wick with 50- 
percent porosity. The binary diffusion coefficient was evaluated from the Chapman- 
Enskog expression (ref. 7) with measured values of the Lennard -Jones parameters for 
the argon -lithium interaction (ref. 8). 

Figure 3 shows the results for startup with a constant heat input Q of 5 kilowatts 
and initial gas pressvires Pj^ of 50, 100, and 200 tor r. The switch on times t^^ for these 
cases are -12.2, -13.0, and -13.9 seconds, respectively. At t = 0 the vapor energy 
transport increases from near zero to slightly less than Q in less than 0. 02 sec- 
ond. In this period of time the vapor pressure p.^(T) increased by only 0. 7 percent, but 
this is sufficient to raise the parameter X (eq. (8)) from 0 to 2. 2. It should be noted that 
it takes about the same length of time to establish a steady -state concentration profile in 
an argon -lithium system of the previously given dimensions at these pressures and tem- 
peratures. Nevertheless, the figures show that in each case startup is initiated at a very 
precise temperature T^ given by equation (36). This fact can be used to eliminate 
startup difficulties in many applications. Since the heat -transfer capability of the heat 
pipe increases with the vapor pressxire p^, the initial gas pressime p^ can be increased 
to the point where startup is not initiated until p.^ is large enovigh to support a heat- 
transfer rate large enough to overcome the thermal load on the condenser. 

Figure 4 shows the results for a startup in which Q is increased in steps to 25 kil- 
owatts with p. = 50 torr. In this way the vapor flow rate at the evaporator exit does not 
reach a very large value until p^ is large enoxigh to sustain it. 
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CONCLUSIONS 


A two-dimensional diffusion analysis has been used to determine the energy trans- 
ported by the vapor between the hot and cold zones of a gas -loaded heat pipe. This 
theory has been combined with a simplified model of the startup process. Numerical re- 
sults from the combined analysis show that very large diffusion rates occur during start- 
up. Such a result cannot be obtained from a one -dimensional axial diffusion analysis 
since the radial flow is an essential part of the process. 

The numerical results also show that startup is not initiated until the vapor pressure 
p.^(T) reaches a precise value which is proportional to the initial gas pressure p^^. 

Hence, startup difficulties can be reduced by selecting a value for p. such that the pipe 
will not start until p^ is large enov^h to support a heat -transfer rate siifficient to over- 
come the thermal load on the pipe. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 11, 1973 
503-25. 
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